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Involute Composite Design Overview 


The strong interaction between material architecture, processing and structural 
performance for nozzle components was described at an earlier NASA symposium for 
laminated involute composites. Since that meeting the Space Shuttle SRM nozzle has test 
fired involute nozzle components and progress has been made in analyzing their 
sensitivity to ply pattern design. The parameters that control ply pattern shape [1,2] also 
control tooling for the manufacture of involute composite structures. In the current 
CAD/CAM idiom these parameters might be called material form features and they 
provide a basis for global composite design sensitivity derivatives. They are not to be 
confused with laminate point design parameters that ignore ply continuity constraints 
present in finite dimension structural components with curvature. We first define the 
involute design problem and illustrate these commonly used approaches for composite 
shell structures. Then analytic sensitivity derivatives are developed and used to analyze 
test rings and cones with maximum stress failure criteria. 
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Involute Exit Cone 
Ply Material Distribution 


The figure below illustrates the orientation of the plies in a conical section of an 
involute exit cone. The intersection of a ply with a plane defined by a constant Z 
coordinate (e.g., curves AD and BC) is an involute curve. Each ply can be mapped 
to an adjacent ply by a rotation of 0 degrees about the axis of symmetry, where 


and N is the number of plies in the involute structure. 
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Involute Design Practice 


Industrial practice in the U.S. for involute ply pattern design at one time was limited to 
either ID or OD start lines partly because of the geometric complexity of the problem. 
Both of these design approaches have one straight edge which also makes layup and 
inspection easier. On the negative side these patterns are more difficult to form to shape 
and do not insure fiber continuity in critical stress regions. An alternative approach used 
by PDA places the start line near the midsurface of the net part to insure fiber continuity. 
Typically a great many patterns are examined by trial and error using CAE tools in 
designing involute composite structures. Low manufacturing risk and high margins of 
safety during a motor firing are the figures of merit. The shape of the component, hence 
weight, is prescribed in most cases and rarely is this shape open to significant change. The 
shape of the ply pattern in contrast is open to wide variations and suggests the need for 
design sensitivity analyses to improve trial and error procedures and ultimately to 
automate the procedure. 



PDA Start-line Ply Pattern 



I.D. Start-line Ply Pattern 
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Involute Design Variables 


The design space includes the ply count and the ply thickness product Nt, the 
helix angle <j> 0 , and up to 6 variables defining the start line, which is a ply 

meridian lying in the r-z plane. These variables determine at each node three 
Euler angles a, y, and <j> which rotate the reference frame into the material 
frame. From these angles a strain transformation matrix is calculated and the 
distribution of these matrices within a finite element is used to calculate the 
element stiffness matrix. This is a global relation independent of finite 
element mesh. 


DESIGN VARIABLES - (X) 


(X) = (Nt, <|>o, R s (Z)) 


DEPENDENT VARIABLES ~ Euler Angles a, y, <j> 


a = a (X) 
y = y (X) 

<1> = 4>(X)J 


[D(a, y, 4>)] 


STIFFNESS MATRIX 


[K] = 


f [B] T [D] T [C m ][D][B] I J I dV 

J V 


995 



Design Sensitivity Formulation 


The direct method of design sensitivity analysis is used. The governing 
equation for linear statics (equation 1) is differentiated to obtain equation 2. 
Equation 3 is obtained by solving equation 2 for {dU/dX}. The remainder of 
the effort is directed toward evaluating {dU/dX}. The Finite element analysis 
already produces the factored stiffness matrix so it is only necessary to 
evaluate the part of equation 3 in parentheses. The element stiffness matrix 
in equation 4 and the element thermal load vector in equation 5 ( only the 
thermal load is sensitive to material geometry orientation ) may be 
differentiated to obtain equation 6. The differentiation is simplified because 
there is no shape sensitivity: the derivative of the strain-displacement 
transformation matrix [B] is zero. The new matrices [Q] and [R] are 
functions of the elasticity matrix [C], the strain transformation matrix [D] 
and its inverse, and the corresponding derivatives, all given in equations 7 
and 8. The integrands in equation 6 are developed in closed form. All 
sensitivity calculations and all finite element analyses are performed by a test 
version of the P/COMPOSITE module in PATRAN. 
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and 

[R] = [C][D] _1 ^ (1D1) 


(8) 
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Vector and Matrix Sensitivity Integration 

Originally the sensitivity integration (equation 6) was performed by 
inserting the {U} vector inside the second integral and making the 
substitution 

{e) = [B]{U} 

to convert the matrix integration to a vector integration. The economy of 
this approach is evident, but finite difference tests show a failure to converge 
manifested by a "plateau" phenomenon for step sizes below a certain 
threshold. The onset of this deviation occurs at a step size that is too large to 
be attributable to round-off error. Because of this error, it was decided that 
matrix integration would be used for all sensitivity calculations. 

The accompanying graph was generated for the helix angle design variable in 
a 439 degree of freedom test cone problem, and the error shown is typical. 


C0MPRR I SQN 0F ANALYTIC AND FINITE DIFFERENCE 
DERIVATIVES WITH RESPECT T0 HELIX ANGLE 
FOR VECTOR AND MATRIX INTEGRATION 
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Formulation of Optimization Problem 


The objective of the optimization is to minimize risk. To this end, the shape 
(and thus the weight) of the part are fixed and the optimization is used to find 
the ply pattern design furthest from the constraint surfaces, subject to 
manufacturing and side constraints. The mathematical formulation of the 
optimization problem is given below. The objective function is a slack 
variable (3 which represents the load margin (i.e., the distance between the 
load index and unity) to be maximized. The slack variable is added to each 
response constraint gj. In addition, there are manufacturing constraints h k 
which do not require the buffer of the slack variable. The Method of 
Feasible Directions algorithm [3] in MICRODOT is used to solve the primal 
form of the optimization problem. Dual methods are not used because the 
number of constraints is much greater than the number of design variables. 
Approximation concepts [4, 5, and 6] are used to formulate the sequence of 
approximate problems. 


(3 — > max 

Subject to gj(X) + (3 < 0 j = l,2,...,m 

h k (X)<0 k = 1,2 £ 

xf^Xi^Xf i =1,2 n 
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Manufacturing Constraints on Design Parameters 

Two of the manufacturing constraints alluded to on the previous page are 
shown in the picture below. The central angle is the angle subtended by a ply 
as it extends from the inner radius to the outer radius, and the arc angle is the 
angle between a ply surface tangent and the circumferential direction. Good 
design practice for nozzle components dictates that the central angle should 
not exceed 120 degrees and the arc angle should not exceed 10-15 degrees. 
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First Optimization Sequence For Involute Ring Problem 

The first optimization problem is an axisymmetric carbon-carbon cylinder 
having 82 degrees of freedom. The design variables are the helix angle <J> and 
the slope m. Five iterations are required to increase the margin of safety 
from .708 at the starting point to 1 .055 at the optimum. 


FIRST 0PTIMIZRTI6N PATH FOR INVOLUTE RING 
INTERNAL PRESSURE LOAD. MAX STRESS CRITERION 
MARGIN OF SAFETY SHONN FOR EACH ITERATION 
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Summary Table for First Optimization Path 

The iteration history of the First optimization path is shown in Table 1 . In the 
course of the optimization the critical stress component varies from fill 
tension in the outer element at the beginning to in-plane shear on the inner 
radius at the optimum. As the optimization proceeds, the constraint tolerance 
is reduced from .03 to the value of .001 required for convergence. The finite 
difference tests on the sensitivity derivatives were used to select move limits 
that would predict the response to within about 10 percent. The optimization 
results indicate that this choice was conservative enough. 
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Second Optimization Sequence For Involute Ring Problem 

A different starting point is used for the same optimization problem. Here 
four iterations are required to increase the margin of safety from .563 at the 
starting point to 1.055 at the optimum. Note that the same optimum is 
reached in both cases. 
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Summary Table for Second Optimization Path 

The iteration history of the second optimization path is shown in Table 2. 
Although the margin of safety shows greater improvement along this path, 
fewer iterations are required. Throughout the optimization the critical 
region is along the inner radius. The critical stress component varies from 
fill tension at the beginning to in-plane shear at the optimum. As before, the 
constraint tolerance is reduced from .03 to the value of .001 required for 
convergence. Again, the move limits appear conservative enough. 
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Involute Test Cone Optimization 


The second optimization problem is an axisymmetric carbon-carbon test 
cone tested and analyzed by Stanton and Kipp [7] and having 439 degrees of 
freedom. The model is subjected to an axial load along the aft rim and is 
constrained axially along the forward rim. The design variables are the 
product of the ply count and the ply thickness Nt, the helix angle <f> 0 , and four 

additional variables controlling the start line. The initial design is the same 
as the final design selected by Stanton and Kipp: it is therefore expected that 
the initial design is close to optimal. Four iterations are, in fact, required to 
increase the margin of safety from 2.653 at the starting point to 2.798 at the 
optimum. 
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Reference and Optimized Ply Pattern 


The reference ply pattern designed by Stanton and Kipp is shown below. The 
warp aligned test cone design with the start line following the midsurface of 
the shell was considered radical when it was first discussed with 
manufacturers. The design goal at that time was to develop the full strength 
of the carbon-carbon material in the critical cone-cylinder transition under 
axial load. That too was considered impractical but both are common 
practice today. 

It can be seen that the shape of the optimized ply pattern in the forward 
region has changed considerably due to changes in the start line. In addition, 
the width of the ply pattern has been uniformly reduced due to the increase in 
the ply count and the ply thickness. 

REFERENCE PLY PATTERN 
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Summary Table For Test Cone Optimization Path 


The optimization history for the test cone problem is given below. 
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Conclusions 


An optimization capability for involute structures has been developed. Its 
key feature is the use of global material geometry variables which are so 
chosen that all combinations of design variables within a set of lower and 
upper bounds correspond to manufacturable designs. A further advantage of 
global variables is that their number does not increase with increasing mesh 
density. The accuracy of the sensitivity derivatives has been verified both 
through finite difference tests and through the successful use of the 
derivatives by an optimizer. 

The state of the art in composite design today is still marked by point design 
algorithms linked together using ad hoc methods not directly related to a 
manufacturing procedure. The global design sensitivity approach presented 
here for involutes can be applied to filament wound shells and other 
composite constructions using material form features peculiar to each 
construction. The present involute optimization technology is being applied 
to the Space Shuttle SRM nozzle boot ring redesigns by PDA Engineering. 


— A design sensitivity capability using global material 
geometry variables has been developed. 

— The number of global variables is insensitive to finite 
element mesh density. 

— The sensitivity derivatives have been used successfully in 
an optimization context. 

— The sensitivity integral calculations in vector form 
yielded a tangible error not shared by the corresponding 
calculations in matrix form. 

— The global variable approach can be applied to other 
composite constructions using material form features 
peculiar to each construction. 
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